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The phase diagram of the random field Ising model on the Bethe lattice with a symmetric dichoto- 
mous random field is closely investigated with respect to the transition between the ferromagnetic 
and paramagnetic regime. Refining arguments of Bleher, Ruiz and Zagrebnov [J. Stat. Phys. 93, 
33 (1998)] an exact upper bound for the existence of a unique paramagnetic phase is found which 
considerably improves the earlier results. Several numerical estimates of transition lines between a 
ferromagnetic and a paramagnetic regime are presented. The obtained results do not coincide with 
a lower bound for the onset of ferromagnetism proposed by Bruinsma [Phys. Rev. B 30, 289 (1984)]. 
If the latter one proves correct this would hint to a region of coexistence of stable ferromagnetic 
phases and a stable paramagnetic phase. 



PACS numbers: 05.45.Df, 05.50.+q, 75.10.Nr, 05.70.Fh 



I. INTRODUCTION 

The random field Ising model (RFIM) has been stud- 
ied extensively in theory [JA as well as in experiment [Q. 
The one-dimensional model ||-[|l6| can be reformulated 
as a random iterated function system (RIFS) for an ef- 
fective field ||, 0> H M. The reformulation leads to an 
iteration of first order whereas standard transfer matrix 
methods lead to iterated function systems of second or- 
der. This considerable simplification allows deep insights 
into the effects of quenched random fields on local ther- 
modynamic quantities. 

Being one-dimensional the Ising chain has no phase 
transitions for finite temperature though. The RFIM on 
the Bethe lattice to the contrary exhibits for not too high 
temperature at least a phase transition from ferromag- 
netic behaviour for small random fields to paramagnetic 
behaviour for large fields J3|, [II |Tfl| . The phase diagram 
is probably even much richer |18|T For T — hysteresis 
effects have been found and investigated in detail |2lJ . 

The Bethe lattice (Cayley tree) is uniquely character- 
ized by the two properties that it is an infinite simple 
graph with constant vertex degree and that it contains 
no loops. It is of order or degree k if the vertex de- 
gree is k + 1. The Bethe lattice of degree k = 1 is the 
one dimensional lattice and the Bethe lattice of degree 
k — 2 the well known binary tree. Because the Bethe 
lattice contains no loops the RFIM on the Bethe lattice 
can be reformulated to a (generalized) RIFS ||, H |T^] 
for the effective field like in the one-dimensional model 
@, §. Therefore, the same powerful techniques as in 
the one-dimensional case can be applied to gain insight 
into the mechanisms driving the phase transition. Nev- 
ertheless the exact transition line in the (T, h) parameter 
plane is still not known. Recently, exact lower bounds 
for the existence of a stable ferromagnetic phase as well 
as exact upper bounds for the existence of a stable para- 
magnetic phase were proved || . We present an improved 
upper bound for the existence of a stable paramagnetic 



phase based on this approach. These bounds are still far 
from the region where the transition is expected though. 
Therefore, we also develop several criteria to detect the 
phase transition line numerically. It turns out that the 
obtained results while being consistent with each other 
disagree significantly with an early result by Bruinsma Q] 
who calculated a lower bound for the onset of ferromag- 
netic behaviour. As Bruinsma's argument rests on the 
differentiability of the density of the invariant measure 
of the RIFS which was only proven for small h and near 
T c there are two possible interpretations. Either Bru- 
insma's bound is not true outside the proven region of 
validity and the transition from ferromagnetic to para- 
magnetic behaviour takes place at the smaller random 
field values found in our numerical results or there is a 
region of coexistence of stable ferromagnetic phases with 
a stable paramagnetic phase implying a phase transition 
of first order in this region. 



The paper is organized as follows. After introducing 
the model and our notations in Sec. II we present the 
improved exact upper bounds for the onset of paramag- 
netism in Sec. III. In section |iy] we give three criteria 
to estimate the transition line between the ferromagnetic 
and paramagnetic regime. The expectation value of the 
local magnetization is calculated directly and we extract 
an estimate for the region of a stable ferromagnetic phase. 
We then study the average contractivity of the RIFS of 
the effective field. This leads to an estimate for the ap- 
pearance of a stable paramagnetic phase for increasing 
random field strength h. The third criterion is the inde- 
pendence of the effective field from boundary conditions. 
It also provides an estimate for the stability region of 
the paramagnetic phase. The implications of our results 
in comparison to Bruinsma's approach are discussed in 
detail in the concluding Sec. |v|. 
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FIG. 1: Bethe lattice of degree k = 2. The solid lines mark 
the part of the lattice denoted by V + and the dashed lines 
the part denoted by V~. The roots of the two subtrees are 
denoted by yo and zo respectively. The thick line shows the 
unique path from a vertex y £ dVs at the boundary to the 
central vertex yo to illustrate the labeling along the path used 
in section ^TJ. 

II. THE MODEL 

The formulation of the RFIM on a Bethe lattice re- 
quires some notations for the underlying lattice. By V 
we denote the set of vertices of the Bethe lattice and 
d(y, z) is the natural metric on the lattice given by the 
length of the unique path connecting y and z. Further- 
more, Vr := {y G V : d(y, yo) < R} denotes the ball of 
radius R around some arbitrarily chosen central vertex 
yo and 8Vr := {y € V : d(y, yo) = R} its boundary, the 
sphere of radius R. In the following it will be useful to 
decompose V into two subtrees V + and V~ with roots 
yo and zo in the way illustrated in Fig. tu. 

Introducing the notation S(y) := {z S 9Vr+i : 
d(y, z) = 1} for the successors of y £ dV R the Hamil- 
tonian of the RFIM on the Bethe lattice reads 

Hn({s v } y ev R ) = Js v Sz ~ E h v s y ~^2 x y s y> 

S/£Vr_i UEVr-i yedVn 

(1) 

where s y denotes the classical spin at vertex y taking 
values ±1, J is the coupling strength, h y is the random 
field at site y and x b y the field at the boundary encoding 
the chosen boundary conditions. We restrict ourselves to 
independent, identically distributed, symmetric dichoto- 
mous random fields, i.e., h y = ±/i with probability ^. 
The canonical partition function 

Z R := eM-(3H R ({s v })), (2) 

{Sy}yeV R 

where /? = {ksT)" 1 is the inverse temperature can be 
reformulated by a method first introduced by Rujan R] 



for the one dimensional RFIM resulting in 

Z R = £exp/3((4«>+A(4?)) %0 + £ B{x^)), 
s yo =±i yev R \{y } 

(3) 

where the effective fields x' z R) are determined by the gen- 
eralized RIFS 

4«>= A & m ) + h y (4) 

with boundary conditions x y a) = x y for y £ dV R . The 
functions A and B are given by 

A{x) = (2/5) _1 ln(cosh/3(x + J)/ cosh/3(a; — J)), (5) 
B{x) = (2/3) -1 ln(4cosh/3(a: + J) cosh/3(ir — J)). (6) 

Note that the upper index (R) of the effective field refers 
to the radius of the sphere where the boundary condi- 
tions are fixed. The partition function in the form (^) is 
a partition function of one spin s yo in two effective fields 
x ( y R) and A(x l z R) ) which are both determined through the 
RIFS (|). The sum in (B) implies that although \A'\ < 1 
for non-zero T, the RIFS is not necessarily contractive in 
contrast to the one-dimensional case. A loss of contrac- 
tivity indicates a phase transition as is explained in more 
detail below. 

Being functions of the random fields h y the effec- 
tive fields are random variables (RVs) on the random 
field probability space and the iteration (|4|) induces a 
Frobenius-Perron or Chapman-Kolmogorov equation for 
their probability measure 

<\X)= \{ H*A#v^)(X-h y ), (7) 

h y =±h z£S(y) 

where Y[* denotes the convolution product of measures, 
X is some measurable set, X — h y := {x — h y \x G X} 
and A# is the induced mapping of A on measures, i.e., 
A#/i(X) :— /i(A^ 1 (X)). The measures of the effective 
fields at the boundary are fixed by boundary conditions, 
e.g. as v { y R) — S x b, the Dirac measure at x b y . Any other 

choice of the RVs x b is also possible though. 

It was proved in [pi that the existence of limiting Gibbs 
measures with finite restrictions compatible with (Q) and 
(^) , cf . |20| , implies the weak convergence of the RVs x y R) , 
i.e., the weak convergence of the measures v { y R) to mea- 
sures v y in the limit R — > oo. For homogeneous boundary 
conditions x h y = x b for all y € V the measures v y are all 
identical and will be denoted by v. 

Before we can present our results on phase transitions 
in the RFIM on the Bethe lattice some more proper- 
ties of the RIFS (Q) and the function A are necessary. 
A{x) is a monotonic function in x. For a given random 
field configuration {h y } yeV + = Wyh} y£V +, <J y = ±, 
we denote the composite function mapping the effec- 
tive fields in dV^ +1 to the effective field at yo by f{ a } R - 
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Here, {ct}r is the tree of k R+1 — 1 symbols ± charac- 
terizing the configuration of the random field and k is 
the degree of the Bethe lattice. These composite func- 
tions are monotonic in the sense that if x h y > x' y for 

all y e dV+ +1 then / Wr ({xJ}) > / Wb ({<}). In the 
same way they are monotonic with respect to the ran- 
dom field, f {a}a ({x b y }) > it ay > °' y for all 

y € V^~. Furthermore, there exists an invariant inter- 
val / = \x*_, x+] with the property that if x y € I for all 
y G dV^ +1 then also f{a} R ({ x y}) € ^ f° r an y random 
field configuration {<r}^. Here, and x* + are the fixed 
points of the composite functions for homogeneous {— } 
and homogeneous {+} configuration of the random field 
respectively. Since A(x) = — A(— x), these fixed points 
are symmetric, x*_ = —x*,. 



III. UPPER BOUNDS FOR THE EXISTENCE 
OF A UNIQUE PARAMAGNETIC PHASE 

In this section we present an exact upper bound for 
the existence of a unique paramagnetic phase in terms of 
the random field strength h. This bound improves earlier 
results in j3j. 

Throughout this section we will use effective fields 
g y := A(xy) in close analogy to the notation in This 
has some advantages in the calculation which will become 
clear below. The iteration (|J) for g y reads 




(for y e dV R ) 
(otherwise) 



(8) 



and we denote the composite functions mapping the ef- 
fective fields {g y } ye9v + i to g yo by f{ a } R . They have the 
same monotonicity properties as the composite functions 

In order to prove the existence of a unique paramag- 
netic phase it is sufficient to show that the RVs g y do 
not depend on the boundary conditions {g y } in the limit 
R — ► oo for any choice of the boundary conditions. We 
use the notation g+ for the effective field at y € V for 
homogeneous boundary conditions g y = g+ in the limit 
R —* oo and g~ for the effective field resulting from 
the corresponding negative boundary conditions g y = g*_ 
where g* + = A(x%) and g*_ = A{x*_). For g+ and g yo 
we use the shorthand notations g + and g~ . Note that 
the dependence of the effective fields on the random field 
configurations is suppressed in this notation. 

Inspired by the proof for the existence of a unique para- 
magnetic phase for the RFIM on the Bethe lattice of 
degree 2 for almost all random field configurations and 
2 < h < 3 in S , we investigate the expectation value 



E{«r}(|fl H 



(9) 



The monotonicity of the composite functions f{ a } R im- 
plies that if this expectation value is zero for the two 



extremal boundary conditions chosen above then it is 
zero for any two sets of boundary conditions. This then 
implies that the RV g yo is independent of the bound- 
ary conditions for almost all random field configurations. 
The goal of this section is therefore to find a criterion for 
the random field strength h which implies that the ex- 
pectation value (^|) is zero. Because of the monotonicity 
of the composite functions f{ a } R w e have g + > g~ and 
thus \g + — g~\ = g + — g~ . Therefore, we consider 



E { „ } (g + -g-)= d V ({<r})(g+ -g-) 



(10) 



E f Mi*}) (9 + i{*})- 9~ ({&})), 
Mb {a} R ={<j} R 



where rj is the product measure of the probability mea- 
sures of the random fields h y = a y h. In the second step 
the integration was split up into a sum of a finite number 
of integrals over sets of configurations with fixed symbols 
{(t}r in Vr and arbitrary {a} S V\Vr. Using the re- 
cursion relation (|J) the integrand can be expressed as a 
function of the effective fields {g^} v edv R on the bound- 
ary of Vr, 

9 + m)-g-({*}) 

= f{5} R -i{{9y}yedV R ) - f{&} R -i({9y}y£dV R ) 



- d 9yf{^} R -A{ 6 z}zeav R }) ■ (gy - g y )• (11) 

y£dV R 

In the second step the mean value theorem has been used 
for f{&} R and S z € [gJ,fl^J are appropriately chosen. The 
partial derivatives in (fTP) are bounded from above by 



nr^A'^ 



, where A'[ R 1 , 



is an upper bound on 



the maximum of A'(x) for x € \q , N , q , N l, the interval 
of possible values of the effective field at the vertic es z; (y) 
along the unique path from y to ya, cf Appendix A 1 for 
details. This bound only depends on {cr} R = and 
hence is independent of the integration. Thus, 



Em(.9 + -.9~) 



R-l 



< 



Hi ax 

Wr yeOVn 1=0 { *} R ={a} R 



n <i 



(12) 

drj({v}){9y - 9y )- 



The remaining integral for each y is bounded from above 
by 2~l y "lE fl where E_r = max. yedVli E {(7 } (g+ - g~), cf 
Appendix A 2. We therefore obtain 



E =E W (<?+-<?-) 



R-l 



< E2 -,^I E Y[A^ y)max ER. (13) 

Mr y&dV R 1=0 

The finite sums commute and as A' {R ! <, is obtained 

z l \V ) max 

with homogeneous boundary conditions the sums X){o-} n 
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FIG. 2: Exact upper bound for the existence of a stable para- 
magnetic phase on the Bethe lattice of degree k — 2 (solid 
line). The bound was obtained as described in the text with 
all random field configurations at R — 4. The dashed line 
is a similar upper bound obtained by considering a sample 
of 10 4 realizations of the random field at R — 11 using the 
complete sum (^). Close to T — the problem is numeri- 
cally unstable; results are presented only for T > 0.1. The 
large dot was obtained for R = 23 using jl5| ) and 10 5 ran- 
dom field configurations. In the shaded region the result of 
j| for the existence of a unique paramagnetic phase applies. 
The grey dashed lines are the ferromagnetic and the antifer- 
romagnetic lines Q], cf. also J3], and the grey dash dotted 
line is Bruinsma's lower bound for the existence of a stable 
ferromagnetic phase (J = 1) 

are identical for all y € 8Vr such that the sum over y 
can be replaced by a factor \8Vr\ = k R yielding 

Eo < KEr, (14) 

where 

K := £ 2-™ & n (15) 

{a} R 1=0 

Because of the translation invariance of the Bethe lat- 
tice these considerations can be applied recursively. This 
implies Eo < K r E r .R. If the factor K is less than 1 
for any parameters (T, h) we immediately obtain Eo = 
-^{cr} (l<? + — 9~\) = as E r .R is uniformly bounded by 
2g* + for all r e N and therefore K r E r . R — ► for r — > oo. 
By translation symmetry this result holds for all g y with 
y E V. As \g + — g~\ > the vanishing expectation even 
implies \g + — g~\ =0 for almost all realizations {a} of 
the random field. 

The reason for using g y instead of x y is now easily 
explained. If we used the effective fields x y instead of g y 
the product over derivatives of A would be from / = 1 
up to R. This gives a less precise estimate because x y 
with y e 8Vr is less restricted than x yo and therefore the 
bound for A'(x y ) with y € OVr is greater than the one 
iorA'(x yo ). 

To apply the criterion obtained above we evaluated K 
on a computer. The calculation time is proportional to 



the number of random field configurations on Vr and 
thus asymptotically grows for, e. g., k = 2 as 2 2 . There- 
fore, the calculation was restricted to R < 4 (for R = 5 
each data point in an array of 20 x 40 points would take 
about 3 days on a Pentium II 350MHz). The solid line 
in Fig. U shows the upper bound for the existence of a 
unique paramagnetic phase obtained for R = A. 

To estimate the results for R > 4 we relied on statis- 
tical methods and sampled random field configurations. 
Doing so it is saving time not to exploit the symmetry 
and use ( |l3[) instead of (|l5|). The resulting bound for 
R = 11 and a sample of lCr random field configurations 
is the dashed line in Fig. g. 

As the obtained bound systematically depends on the 
radius R it is tempting to try to extrapolate to R = oo. 
For T — 1.2 we obtained an extremely good fit for the 
data sampled at R = 5, . . . , 23 using h c (R) = a + bR c . 
However, the result was a = h c (oo) « which is not 
realistic. Other fits, e. g. omitting data for small R or 
using different functional forms, yielded values between 
—0.9 and 0.45. We suspect that a naive choice of the 
functional form of h c (R) does not allow realistic extrap- 
olation results for these bounds in the case of the Bethe 
lattice. 



IV. ESTIMATES FOR THE TRANSITION LINE 
A. Direct calculation of the magnetization 

Even though the bounds presented in the preceding 
section considerably improve former analytical results 
they are still far away from the region where the phase 
transition from paramagnetic to ferromagnetic behaviour 
is suspected. In [Q Bruinsma claimed to have found a 
lower bound in h for the existence of a ferromagnetic 
phase which is in the relevant parameter region, cf. Fig. |^. 
To check this bound and to get a good numerical approxi- 
mation of the transition line we developed several numer- 
ical criteria for the existence of a ferromagnetic phase or 
the existence of a stable paramagnetic phase. 

The most obvious criterion for the existence of a fer- 
romagnetic phase is a non-vanishing expectation value 
for the magnetization for small but non-zero boundary 
conditions. The expectation value for the local magneti- 
zation at the spin in the center is given by 

m := E w (s yo ) 

- J du{x) dv{y) tanh (0(x + g(y))) (16) 

where (•) denotes the thermodynamic average, E/^j the 
expectation value with respect to all random field con- 
figurations and v is the limit measure of the effective 
field for homogeneous boundary conditions x y = x b for 
all y S V in the limit R — > oo. To approximate v we 
generated a large number of random field configurations 
on a finite region Vr and calculated the corresponding 
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FIG. 3: Estimates for the transition line between between 
para- and ferromagnetism. The dashed line separates regions 
where the average magnetization decreases or increases with 
the distance from the boundary for homogeneous boundary 
conditions x b = 0.01 obtained by comparing the configura- 
tional averages over 4 ■ 10 up to 64 ■ 10 5 samples at distances 
9 and 13 from the boundary. The dotted line gives the bound- 
ary of the region in which the iteration (^) is contractive on 
the average above and non-contractive below. It was obtained 
by sampling 10 5 field configurations and evaluating ([l8]) with 
R — 13 and Ri = 4. The solid line was obtained using (p^) 
at R = 4. The big dots were obtained by evaluating (BJ) for 
R = 20 and between 10 5 and 2.4 • 10 6 random field config- 
urations. The grey lines are as in Fig. |^. The grey shaded 
region marks the region all our numerical results (including 
considerably more than shown explicitly here) are contained 
in. (k = 2, J = 1) 



B. Average contractivity of the RIFS 

For zero boundary conditions there is a paramagnetic 
state for any temperature T and random field strength 
h. The stability of this state is tied to the average con- 
tractivity of the RIFS (||). If it is globally contracting 
the paramagnetic state is stable and unique. If it is at 
least contracting on the average for some interval around 
zero, the paramagnetic phase is stable but the existence 
of other stable phases is not a priori excluded. The in- 
vestigation of the contractivity of the RIFS was first pro- 
posed in [l7| . 

We generated a set Hr of random field configurations 
{cj-R on a finite ball Vr and calculated the image of 
a small initial interval I b = [— x b ,x b ]. Because of the 
monotonicity of f{ a } R the image of this interval at vertex 
y is I y = [x y R ^(— x b ), x y R) (x b )]. To estimate the average 
contractivity of the RIFS we compared the average length 
l/k Rl J2y£dv R \Iy \ at tne vertices y £ 8Vr 1 to the length 
\I yo \ at the central vertex yo. As the effective fields at 
all y £ 8Vr x contribute to the effective field at yo we 
consider the average interval lengths at vertices in 8Vr 1 
instead of individual values. To minimize the influence of 
the somewhat arbitrary choice of the initial interval the 
comparison was performed for R\ <C R. 

There are two ways of performing the comparison. Ei- 
ther one first averages over the lengths \I y \ at all y £ 0Vr 1 
then calculates the quotient of \I Vo \ and this average 
length in 8Vr 1 and average over the sample Hr of ran- 
dom field configurations at the end, 



effective field x y R) . The obtained values were then sorted 
into small boxes of length e. The resulting histogram was 
used as an approximation of v yo (bi) —: where hi is the 
i-th box. Explicitly this yields for the magnetization 



ViV i ta - uh P(x i + g(y.j)). 



(17) 



where the points Xi and yj were chosen as the center of 
box i and j respectively. 

Assuming that the magnetization in the center varies 
monotonically with the radius R of the finite volume Vr 
one would expect to observe a monotonically increasing 
magnetization in the ferromagnetic regime and a mono- 
tonically decreasing magnetization in the paramagnetic 
regime for increasing R. Therefore, the dashed contour 
in Fig. |^ which divides the regions in which the numerical 
estimate of the magnetization is increasing or decreasing 
with increasing R is a good estimate for the transition 
line. This type of estimates only slightly depends on the 
chosen boundary condition and iteration depth but the 
obtained results significantly disagree with Bruinsma's 
bound. 



yedV Rl \ 2 V\ 



(18) 



Or one first averages \I y \ over all y £ 8Vr 1 and all ran- 
dom field configurations as well as \I Vo \ over the same 
random field configurations and calculates the quotient 
at the end, 



(19) 



The two averaging procedures ( pL8| ) and (19) yield identi- 
cal results and thus obviously are equivalent. 

If the images of the initial interval contract on the 
average for a finite iteration of (|^) we expect complete 
contraction to length zero for infinite iteration. This cor- 
responds to a stable paramagnetic phase. Therefore, the 
contour in the (T, h) parameter plane at which the aver- 
age quotient of band lengths switches from greater than 1 
below to less than 1 above is an estimate for the stability 
region of the paramagnetic phase. The resulting esti- 
mated transition line is shown for R = 13, Ri = 4, the 
initial interval I b = [-0.01,0.01] and 10 5 random field 
configurations as the dotted line in Fig. ||. Again, the 
agreement of the obtained results for various boundary 
conditions and iteration depths is satisfactory but there 
is a large deviation from Bruinsma's line. 
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Independence of the effective fields from 
boundary conditions 



1. 



+ - + - + - 
+ - + 



2. 



+ + - + 
+ 



A related criterion for the existence of a stable para- 
magnetic phase is the independence of t he effective field 
from boundary conditions. As in Sec. Ill we use the effec- 



tive fields g y R) rather than x y R) . We consider boundary 
conditions {g y }y£dv R taking values in a small interval 
[— g b ,g b ]. Through the iteration with (||) the effective 
fields g y R> are functions of the boundary conditions 

9? = fr.H-nMvM}) (20) 

where the function f{ a } R _ n „ 1 (y) has k R ~ n arguments for 
y S dV n and it is the identity if R = n. For simplicity 
and without loss of generality we restrict the following 
discussion to g y R) . The boundary conditions can be writ- 



To 



ten as g b y — g b ' g b where g y takes values in [—1,1]. 
investigate how the effective fields depends on the bound- 
ary conditions we consider the expectation value of the 
derivative of g y R) with respect to g b , the strength of the 
applied boundary condition 



E 



dg b 



y&dVR 

< E ^u-M A 'KU xh )) ( 21 ) 

v eav R i=o 

where x (R> ^(x b ) denotes the effective fields along the 
unique path from y to yo with homogeneous boundary 
conditions x b = A~ 1 (ff fc ). Now one can estimate similar 



as in Appendix A 1 



=: A 



/(«) 
z l(y)i 



(22) 



As the boundary conditions are homogeneous this implies 



E 



dg 



b y yo 



<k J 



R-l 
ill* 



(R) 

2 ife)max 



(x b ). (23) 



If the right hand side vanishes for R — > oo the effective 
field g yo is on the average independent from boundary 
conditions taking values in [— g b ,g b ]. By determining the 
parameter region in which the right hand side of (|23| ) 
vanishes for R — > oo we therefore get an upper bound on 
the emergence of a stable paramagnetic phase. As our 
calculations are limited to finite R, convergence to zero 
is assumed if the obtained values of (j||) for R > are 
smaller than 1 which is the value for R = 0. 

For the Bethe lattice of degree k = 2, radius R = 4 and 
g b = 0.01 the right hand side of (p3|) was evaluated. The 



3. 



+ - + - 
+ 



4. 



+ - + 
+ 

+ 



FIG. 4: Four equivalent chessboard configurations at R — 
4. The second configuration is obtained from the first by 
permutation of two subtrees of a vertex in the sphere dV2, 
the third one by permutation of two subtrees of a vertex in 
the sphere dVi and the last one by permutation of the two 
subtrees of the root itself. 



contour between values smaller than 1 above and greater 
than 1 below is shown as the solid line in Fig. ^. For 
R > 4 we again relied on sampling random field configu- 
rations instead. The resulting transition lines for various 
iteration depths and boundary conditions are compara- 
ble to the results of the preceding two sections and are 
all contained in the grey region in Fig. [| 

If we consider the derivative of the effective field at y$ 
in the case of boundary conditions g = we have 



E 



HI 



R-l 



dg 



b ({0})=E Wa _ l e U A '^ 

yedVn 1=0 



.(«) 

z i{y) 



({0})). 
(24) 



If this derivative does not tend to zero for some param- 
eters (T, h) and R — > oo there is no stable paramag- 
netic phase. By determination of the parameter region 
in which this is the case we get a lower bound on the 
emergence of a stable paramagnetic phase. The numeri- 
cal results are the large dots in Fig. ||. 

V. DISCUSSION 

In order to interpret the discrepancies between Bru- 
insma's result and our numerical investigations we briefly 
review Bruinsma's argument (J] in our language. The 
probability measures v y of the effective fields x y are fixed 
points of the Frobenius-Perron equation, Eq. (Q). They 
can be approximated by finite iterations of some ini- 
tial probability densities (boundary conditions) v b for 
y G 8Vr. If the support of v b is a subset of the in- 
variant interval /, the support of u Vo is a subset of the 
images of I by functions f{ a } R - These images are called 
bands. The left and right boundary of the bands are the 
effective fields corresponding to homogeneous boundary 
conditions x b = x* + and x b = x*_ for y S 8Vr, respec- 
tively. The investigation of the structure of the set of 
bands has proved to be a powerful tool in the treat- 
ment of the one-dimensional random field Ising model 
g §, 0, H [H [14[ H 0. In contrast to the one- 
dimensional case the bands are highly degenerate here, 
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FIG. 5: Approximations of the invariant measure density p 
obtained from four-fold application of (Q) to the equipartition 
p b — ■ 1/ on the invariant interval /. There are no max- 
ima at ±h in a) whereas the two maxima at ±/i in b) are 
already so pronounced that differentiability of p is question- 
able. The random field strength in b) was chosen such that 
the point (T, h) is very close to Bruinsma's bound, (k = 2, 
= J=l) 



i.e., different configurations of the random field result 
in the same band. This is due to the invariance of the 
model with respect to permutations of subtrees for ho- 
mogeneous boundary conditions. The most degenerate 
bands correspond to the two chess-board configurations, 
cf. Fig. ^, of the random field with +h or —h at yo re- 
spectively. There are 2 2 _1 equivalent random field 
configurations in the case of the Bethe lattice of degree 
k = 2 and radius R. As the total number of configura- 
tions is N = 2 2 _1 the most degenerate bands have a 
weight of 2 2R ~ 1 - 1 /2 2R - 1 = 2- 2R ~ 1 ~ l/y/N. The bands 
with the least weight are the bands corresponding to ho- 
mogeneous +h or — h random field configurations. They 
have the weight 1/N. The weights of all other bands are 
distributed between these values. 

Bruinsma used boundary conditions v~ = 5 x b for some 
x b e K and iterated (f?t). He only considered the lowest 
and highest weight terms corresponding to the least and 
the most degenerate bands. The highest weight terms 
obey a recursion relation. The fixed points of this relation 
can be calculated and it is straightforward to determine 
for which temperatures T and random field strengths h 
they are symmetric to the origin. Proving differentiabil- 
ity of the density p of the invariant measure v in a neigh- 
bourhood of T — T c and h = 0, Bruinsma concluded that 
an asymmetric position corresponds to asymmetric max- 
ima of p of non-zero weight and therefore to the existence 
of a ferromagnetic phase. 

The symmetric position corresponds to complete con- 
traction of the most degenerate bands such that the 
asymmetric boundary condition has no effect in the limit 
of infinite iteration. The asymmetric position on the 
other hand occurs if the most degenerate bands do not 
completely contract. Seen this way the argument above 
is o ur cri terion of average contractivity of the RIFS in 
Sec. IV B except that it is restricted to the contractivity 



local maxima in a differentiable measure density really 
is given by the most degenerate bands. For small h this 
actually seems not to be the case, cf. Fig. ||a. As the 
maxima are at ±/i and therefore close to zero for small 
h it is difficult to argue based on numerical data though. 
The example in Fig. |^b shows however that the maxima 
are present for sufficiently large h. 



Secondly, the differentiability of the invariant mea- 
sure density has been proved only in a neighbourhood 
of T — T c and h = whereas for large h or small T, 
the measure density p is clearly not differentiable. It is 
unclear whether it is differentiable in the region of the 
lower bound, cf. Fig. @b. 



The disagreement of our numerical work with Bruins- 
ma's lower bound therefore allows two interpretations. 
Either Bruinsma's bound is not true outside the proven 
region of validity because the most degenerate bands are 
not a sufficient indicator for the symmetry of p when 
the measure density is not differentiable. Or in the re- 
gion between our upper bounds for the existence of a 
stable paramagnetic phase and Bruinsma's lower bound 
for the onset of ferromagnetism a stable paramagnetic 
phase coexists with the — also stable — ferromagnetic 
phases. This would imply the existence of a first order 
phase transition and hysteresis loops depending on the 
strength of the random field in contrast to the hysteresis 
at T = |2lj] which depends on the homogeneous offset 
of the random field. 



In this paper we improved exact upper bounds for 
the existence of a unique paramagnetic phase in Sec. Ill 
which is a further step towards the exact determination of 
the phase diagram of the RFIM on the Bethe lattice. Fur- 
thermore, we presented numerical work leading to various 
estimates for the actual phase transition line. The direct 
calculation of the expe ctation value of the local mag- 
netization in Sec. IV A , the investigation of the average 
cont ractiv ity of the RIFS (|4|) at large iteration depths in 
Sec. IV B and the numerical calculation of the derivative 



of one specific band instead of the average contractivity. 

There are two problematic points in the reasoning 
above. Firstly, it is not clear whether the location of 



of the effective field with respe ct to the strength of the 
boundary condition in section IV C provided estimates 
for the stability region of the paramagnetic phase. All 
results are in good agreement while all disagreeing with 
the earlier result of Bruinsma Q . This disagreement mo- 
tivates further investigations whether the bound for the 
onset of ferromagnetism given in Q needs to be recon- 
sidered or whether there really is a coexistence region for 
stable ferromagnetic and a stable paramagnetic phase. 



The work was partially supported by the Cusanuswerk 
and the DFG (Graduiertenkolleg "Quantenfeldtheorie"). 
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APPENDIX A 



1. Bound on the partial derivatives in (11) 



The partial derivatives in ( |11[) are given by 
9g y f{a} R ^ 1 ({Sz} zedv R }) 

R-l 

= n A '(f{v}R-i-i(zi(y))({ £ z}zedV R )), (Al) 

where e z = A~ 1 (S Z ), zi(y) £ dVi are the vertices along 
the unique path from y to y$, cf. also Fig. |[ and 
{a} R-i-i(zi(y)) are the signs of the random field con- 
figuration on the subtree of depth R — 1 — I with root 
zi(y). The terms f{B} R . 1 . l {z l ( v )){{ s z}zedV H ) are effec- 
corresponding to boundary conditions 



tive fields 



r (H) 

'My) 



{x b z = e z }z£dv R - We write x^]Ae) for these fields and 
x zf( y \( x+ ) an d x ^(y)( x ~) ^ 0r ^ ne corresponding effective 
fields with boundary conditions x z — xf and x\ = x~ 
for z <G 8Vr. We then can estimate 



<max{A'«> Je')) 



max 
min 



Z i J i 

•'' \ X zi(y)( X ^ X z,(y) 



{A'(x) 

A'(min{x^ y) (x+),0})} 



x z 

.(B) 



G dV R ] 
(* + )]} 



(A2) 



In the last step we used that the maximum of A' in an 
interval [a, b] is at a if a > 0, at b if 6 < and at zero 
in all other cases. As the effective fields can never be 
larger than x* + and never smaller than x*_ we can for 



z e 8Vr estimate x z > h z + kA(x*_) 



^{x*_) and 



xf < h z + kA(xl) = x z n+1) (x* + ). This allows to replace 



x + and x in the argument of x^] y \ i n ( A2) and with 

x( z?(y)K +1) ( x *±)) = < ( ;>±) we get 



A'(x z *l y) ( s )) < u^{A\m^{x^( x -)m 
A'(mxn{x^(x* + ),0})} 



A' 



(A3) 



Inserting ( |A3j ) into (Al) then yields 

R-l 

d 9y h} R ^z}zeov R })<l[Kly)r 



(A4) 



2. Bounds on the integrals in (12) 



Using the independence of the RVs g y of the signs 
{&z} z <EV R \{y} an d denoting the number of vertices in Vr 
by \V R \, i.e., \Vr\ = (k R+1 - l)/(fc - 1), one obtains 

dv({o-})(9y ~9y) 

{a} R = {a} R 

= 2- |y « l+1 J dr 1 ({a})(g+-g;) 

Oy=<Jy 

= 2- 1 v * I E (g+ - g- \a y =a y ). (A5) 

The function A is antisymmetric which implies 
<7+ ({— a}) = — g~ ({a}) and therefore 

E(.9+ - g; \<r v = +) = E(.9+ - 9y I ^ = -), (A6) 
implying 



±E(.g+ - g~ | <7„ 



= 5 E (.9,} - 5jT I °V = +/ > 2 
= E(.g+ - | a, = &), (A7) 

for any cr 6 { — ,+}. Setting 

E fi := max E M ( 5 + - #~ | a y = <r y ) 

y<EOV R 

= m ax E {ff} 0+ -3"), (A8) 
this finally yields 

d7 ? ({^})( 5 +-p-)<2-l^lE ii . (A9) 



{*}ii={o-}i? 
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